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Measurement of entanglement remains an important problem for quantum information. We 
present the design and simulation of an experimental method for entanglement estimation for a 
general multiqubit state. The system can be in a pure or a mixed state, and it need not be "close" 
to any particular state. Our method, based on dynamic learning, does not require prior state recon- 
struction or lengthy optimization. Results for three-qubit systems compare favorably with known 
entanglement measures. The method is then extended to four- and five-qubit systems, with relative 
ease. As the size of the system grows the amount of training necessary diminishes, raising hopes for 
applicability to large computational systems. 
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I. INTRODUCTION 

Because entanglement is crucial to most quantum in- 
formation applications, convenient and accurate methods 
for its measurement and calculation are essential. Alas, 
though considerable work by a large number of very intel- 
ligent people has made considerable progress (see, e.g.,[T]- 
IT2"]). we have as yet no general method, of universal ap- 
plicability, much less one which is easy to implement. 

For 2-qubit systems the problem is fairly well under- 
stood, from more than one point of view. For exam- 
ple, the "entanglement of formation" [TJ [5] approach is 
based on the idea of construction: how much entangle- 
ment is necessary in order to reconstruct this particular 
state? This gives us a general measure as an analytic 
formula, so that, if one knows the density matrix, one 
can easily calculate the entanglement. However, gener- 
alizations to larger systems are difficult [5]. Various geo- 
metric methods [H|S] use an extremely appealing idea of 
"closeness" to the subspace of separable states: one con- 
structs some distance measure in the state space, then, 
performs a minimization calculation to find the distance 
of closest approach. Product states, of course, have zero 
distance to that space, and have therefore zero entangle- 
ment. Many entanglement witnesses[6] are based on this 
same idea. The drawback is, of course, the optimization 
calculation, which rapidly becomes impractically difficult 
for large systems. In addition, these approaches assume 
that one knows the density matrix. One might not. And 
reconstruction of the density matrix, like optimization, is 
not easy to generalize to large systems. Most witnesses 
also require that the state of the system be "close" , in 
some sense, to a given, known state[7]. 

Clearly it would be of great interest to have a general, 
experimental method for estimation of the entanglement, 
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for which one would not have to know or to reconstruct 
the density matrix, which would not require lengthy op- 
timization procedures, and which could be applied to any 
state. We present here such a method. We demonstrate 
that our method works for large classes of states of 2- 
qubit, 3-qubit, 4-qubit, and 5-qubit systems, and show 
that extension to larger systems is clearly achievable. 

II. DYNAMIC LEARNING: QUANTUM 
NEURAL NETWORK (QNN) 

We begin with the Schrodinger equation: 

t - &™ w 

where p is the density matrix and H is the Hamiltonian, 
whose formal solution[I3] is 

p(t)=exp(tLt)p(0), (2) 

where L is the Liouville operator. We consider an N- 
qubit system whose Hamiltonian is: 

N N 
H = 22 K a (J xa + £a&za + ^ CapCza&zP (3) 
Q=l a^P=l 

where {a} are the Pauli operators corresponding to each 
of the qubits, {K} are the tunneling amplitudes, {e} are 
the biases, and {£}, the qubit-qubit couplings. We choose 
the usual "charge basis ", in which each qubit's state 
is given as or 1; for a system of N qubits there are 
2^ states, each labelled by a bit string each of whose 
numbers corresponds to the state of each qubit, in order. 
The amplitude for each qubit to tunnel to its opposing 
state (i.e., switch between the and 1 states) is its K 
value; each qubit has an external bias represented by its 
e value; and each qubit is coupled to each of the other 
qubits, with a strength represented by the appropriate 
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£ value. Note that, for example, the operator a X A = 
a x (£) I ... <g> /, where there are (N-l) outer products, acts 
nontrivially only on qubit A. 

The parameters {K, £,£} direct the time evolution of 
the system in the sense that, if one or more of them is 
changed, the way a given state will evolve in time will 
also change, because of Eqs. 1-3. This is the basis for 
using our quantum system as a neural network. There 
is a mathematical isomorphism between Eq. [2] and the 
equation for information propagation in a neural net- 
work: <t>out V ut = F W (j)i n put, where 4> outp ut is the output 
vector of the network, §i npu t the input vector, and i*V 
the network operator, which depends on the neuron con- 
nectivity weight matrix W. Here the role of the input 
vector is the initial density matrix p(0), the role of the 
output is the density matrix at the hnal time, p(tf), and 
the role of the "weights" of the network is played by the 
parameters of the Hamiltonian, {K, £,C}> ah of which 
can be adjusted experimentally as functions of time|14j. 
By adjusting the parameters using a neural network type 
learning algorithm we can train the system to evolve in 
time to a set of particular final states at the final time tf, 
in response to a corresponding (one-to-one) set of given 
inputs. Because the time evolution is quantum mechan- 
ical (and, we assume, coherent), a quantum mechani- 
cal function, like an entanglement witness of the initial 
state, can be mapped to an observable of the system's 
final state, a measurement made at the final time tf. 
Complete details, including a derivation of the quantum 
dynamic learning paradigm using backpropagation[15 in 
time [To], are given in [17] , 

The time evolution of the quantum system is calculated 
by integrating the Schrodinger equation numerically in 
MATLAB Simulink, using ODE4 (Runge-Kutta) , with a 
fixed integration step size of 0.05 ns[18]. The system was 
initialized (prepared in) each input state, in turn, then al- 
lowed to evolve for 300 ns. A measurement is then made 
at the final time. All of the parameters {K, e, (} were 
taken to be functions of time; this was done in the simu- 
lation by allowing them to change to a different constant 
value every 75 ns (i.e., four "time chunks".) Discretiza- 
tion error for the numerical integration was checked by 
redoing the calculations with a timestep of a tenth the 
size; results were not affected. For the backpropagation 
learning, the output error needs to be back-propagated 
through time [16], so the integration has to be carried out 
from tf to 0. To implement this in MATLAB Simulink, 
a change of variable is made by letting t' = tf — t, and 
running this simulation forward in t' in Simulink. 



III. WHY THIS IS NOT QUANTUM CONTROL 

More generally, a nonlinear system can be written as 
x = f(x,u,p) (4) 

y = g(x, u) (5) 



where x is the state, f is a nonlinear function, g is an 
output measure, u is the system external control input, 
and y is the system output, x, u, and y are all functions 
of time. Notice from Eq. [T] that our quantum system 
is of the form x = Ax, and the (complex) operator A 
contains the parameters p = {K,e,C}. In general p can 
be a set of a priori fixed system parameters, (e.g. the 
mass, spring constant, and damping coefficient, in a mass 
spring damper system.) u is a control input that is not 
fixed ahead of time, and can be calculated on the fly and 
varied as the system is running, i.e., for every different 
run of the system, u can be recalculated based, usually, 
on the current state of the system x(t) as it is going from 
x(0) to x(tf). So u can be different for different initial 
conditions (x(0) values). While the set of parameters p 
can be functions of time as the system runs, they are 
fixed at the same values every time the system is run 
going from any x(0) to x(i/), once determined via the 
network learning process. 

A control problem, then, would be to calculate u on 
the fly over a period of time in order to make the system 
states x or output y behave in a certain way. An optimal 
control problem would be to choose u to minimize some 
cost function J(x,u) or J(y,u). Given initial conditions 
x(0) and u, the system can be integrated from t=0 to 
t=tf so that 

x(t f ) = F{x(0),u,p) (6) 
y(tf)=g(x(tf),u(tf)) (7) 

Our quantum system, on the other hand, has no external 
input, so it is written 

x(t f ) = Q(x(0), P ) (8) 
y(t f ) = m(x(tf)) (9) 

where now Q defines the quantum system and m is some 
output measurement made on the system at the final 
time. The quantum computation problem is to find val- 
ues for the specific fixed parameters p such that the quan- 
tum system approximates a given function s = h(f) at 
a finite number of points, m, in the domain, that is, 
y{tf) = §i when x(0) = fj i = l,m. 

This is not a control problem, but a problem of find- 
ing a finite dimensional (with dimension equal to the 
number of parameters in p) subspace approximation of 
the function h. Interpreted as a quantum neural net- 
work, the weights p are found through a learning algo- 
rithm, where [rj,Si] i=l,m are the training pairs in a 
training set, and the cost function to be minimized is 
J = [§i — yi(tf)] 2 for each and every i=l,m, where from 
Eqs. 8-9, yi(t/)=m(Q(xi(0),p)) is the output for the sys- 
tem started in initial condition Xi(0)=ii. 

IV. ENTANGLEMENT 

Almost all the information contained within a system 
of 3 qubits is contained within the 2-particle reduced 
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TABLE I. Input matrices for training the pairwise quantum 
neural network (QNN) entanglement witnesses. Each column 
is an input state, showing relative amounts of each of the basis 
states. The N-qubit system was trained on this set of four 
( N\ 

for each of the I ^ I pairs; this constituted the "pairwise" 
training. 



Input 


Bell 


Flat 


Corr. 


P 


|00) 


1 


1 





1 


|01) 





1 





1 


|io> 





1 


0.5 


1 


111) 


1 


1 


1 






density matrices 8J. Thus it makes sense to train our 
entanglement measure, to begin with, to the set of pair- 
wise entanglements. Fortunately in previous work[T7] 
we have already found a set of parameters that success- 
fully maps the input state of a two-qubit system to a 
good approximation of the entanglement of formation, 
via the qubit-qubit correlation function at the final time, 
(a z A(t f)czB(t f)) 2 . The parameters were found by train- 
ing with a set of just four inputs, as shown in Table [TJ a 
fully entangled state, a "flat" state (equal amounts of all 
basis states), a correlated product state, and a partially 
entangled state. 

For the three-qubit state ABC there are three pos- 
sible pairs: AB, AC, and BC; and therefore three dif- 
ferent pairwise entanglements to train. We therefore 
used|19| a set of twelve training pairs: the four in Ta- 
ble IT] for each pair, outerproducted with |0) for the 
unpaired qubit. So, for example, the BellAc state 
was taken to be ^(|000) + 1 10 1) ) , and was trained to 
give an output of 1 for the network output measure 
Oac = (&zA{tf)<T z c(tf)) 2 and zero for each of the out- 
put measures Oab = (&zA(tf)o'zB(tf)} 2 and Obc = 
(&zB(tf)&zc(tf)) 2 - Table |ll| shows the target outputs 
and average calculated outputs after training; detailed 
procedures and results are in |19j . Note that we report 
here the averages; e.g., the Bell a p column numbers are 
the averages over a/3 = AB, AC, and BC. ("Non-a/3" 
refers to the other two outputs for each pair; there are 
six of these numbers.) 

After the pairwise training was accomplished for the 
3-qubit system, we then trained the residual[9] entangle- 
ment, as well. We did this by expanding the training set: 
adding one additional training pair, the three-way en- 
tangled Greenberger-Horne-Zeilinger (GHZ) [3] state for 
three qubits, ^(|000) + |111)), for which we specified tar- 
gets of zero for all three pairwise correlation functions, 
and 1 for the three-point correlation function O abc — 
(&zA(tf)<T Z B(tf)crzc(tf)) 2 - The bottom row in Table [n] 
displays these targets and calculated average outputs. 

Once trained, the parameters that are found can be 
used to evaluate the entanglement of any desired state. 
(This will tell us whether we actually have an estimate 



of the entanglement, or have just been curvefitting the 
proverbial elephant.) Testing was therefore done on a 
large number of states not represented in the training 
set, including fully entangled states, partially entangled 
states, product (unentangled) states, and also mixed 
states. Detailed results were presented in [19]; while er- 
rors were significantly larger, it was clear from the output 
matrix that there is definite separation, in two senses: 
first, that it is easy to see where the pairwise entangle- 
ment is (e.g., to distinguish between a state with AB 
entanglement and one with BC entanglement); and, sec- 
ond, that it is easy to differentiate among unentangled, 
partially entangled, and fully entangled states. In other 
words, while this is not, strictly speaking, an entangle- 
ment measure, it is quite a good witness. 

It was natural to take this one step further, to a sys- 
tem of four qubits. Again, we started from the trained 
parameters for the smaller (3-qubit) system, and simply 
copied the new sets necessary (one tunneling parameter, 
one bias parameter, and three coupling parameters) from 
the trained sets already found. Results are shown in Ta- 
bles [nlj |IV[ and [V] Again, we worked in stages: start- 
ing with a training set of 24 for the pairwise training (6 
pairwise Bell states, 6 pairwise flat states, 6 pairwise cor- 
related unentangled states, and 6 pairwise partly entan- 
gled). Once those were trained, we added four more pairs 
corresponding to the four distinct 3-way GHZ states; af- 
ter those were trained, we added also the single 4-way 
GHZ state. Testing was then done, as with the 3-qubit 
system, on a similar (but expanded) set of states. Errors 
were slightly larger but the results maintained good sep- 
aration. Buoyed by confidence in our method we then 
successfully expanded to a 5-qubit system, using an ex- 
actly similar procedure. Training results are shown in 
Tables [VT VII and VIII Final values for the parameters 
are in TableTm 



V. RESULTS AND DISCUSSION 

Training of the three-qubit system required consider- 
able time: 5000 epochs (passes through the entire train- 
ing set) to get pairwise entanglement from the trained 
results for the two-qubit system [19]. But training of the 
four-qubit system was (unexpectedly) easy. The training 
sets were, of course, larger: there are only three ways 
to entangle a pair of qubits, and only one way to en- 
tangle a triple, for the 3-qubit system, while the number 
of ways for the 4-qubit system are 6 and 4, respectively, 
with one quad (and so on, following the binomial coeffi- 
cients.) However the 4-qubit system required much less 
training: only 100 epochs at each stage. The 5-qubit sys- 
tem was trained for 100 epochs at each stage, but 98% 
of the change in the parameters, as well as in the error 
reduction, had been accomplished by less than half that 
time, as shown in Figure [l] Successively larger systems 
still need to be trained to decrease initial errors; yet the 
amount of training necessary seems to diminish with sue- 
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TABLE II. Stages 2 & 3: Three-qubit QNN entanglement targets and calculated average outputs. Each column is an input 
state corresponding to the respective column in Table]!] each row corresponds to a target function for the specified entanglement 
function. Stage 2 started from the trained 2-qubit system's parameters (stage 1) [TT] , and trained the three-qubit system to the 
set of pairwise entanglement training pairs(12 training pairs, 3 outputs per pair); Stage 3 used the pairwise set plus the GHZ 
state (13 training pairs, 4 outputs per pair). The double vertical line separates the states included in the training set from 
those only tested on. Nonzero target-training pairs are boldfaced for easy comparison. RMS errors were: 1.2 x 1CT 3 for stage 
2, and 1.8 x 1CT 3 for stage 3. 
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0.0014 
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0.0081 
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0.0037 


0.5735 



Targets: 


a/3 


1 








0.44 
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Trained 
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0.9948 


0.0001 


0.0009 


0.4389 






non-a/3 


0.0005 


0.0012 


0.0004 


0.0022 


0.0035 




ABC 


0.0013 


0.0006 


0.0001 


0.0012 


0.9883 



TABLE III. Stage 4: Four-qubit QNN entanglement targets and calculated average outputs. Stage 4 trained the four-qubit 
system to the set of 24 pairwise entanglement training pairs. RMS error = 2.2 x 10~ 3 (24 training pairs, 6 outputs per pair). 
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TABLE IV. Stage 5: four-qubit Q NN entanglement targets and calculated average outputs. Stage 5 added the four 3- way GHZ 
states to the training set of Table III RMS error = 1.9 x 10 -3 (28 training pairs, 10 outputs per pair). 
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0.0014 
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0.0014 


0.0012 


0.0000 


0.0036 


0.0010 


0.7649 



cessive training stages. 

Figures [2] and [3] make this point graphically. They 
show the evolution of the tunneling, bias, and coupling 
parameters, from initial training of the 2-qubit system, 



through each stage of the training of the 3-qubit, 4-qubit, 
and 5-qubit systems. Initially the QNN's learning entan- 
glement requires large changes in the parameters, but 
as the size and complexity of the system is increased, 



5 



TABLE V. Stage 6: four-qubit QNN entanglement targets and calculated average outputs. Stage 6 added the one 4-way GHZ 
state to the training set of stage 5. RMS error = 2.2 x 10 -3 (29 training pairs, 11 outputs per pair). 
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0.0014 
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TABLE VI. Stages 7 & 8: five-qubit QNN entanglement targets and calculated average outputs. Stage 7 trained the five-qubit 
system to the set of 40 pairwise entanglement states for the 5-qubit system; Stage 8 added the ten 3-way GHZ states. RMS 
errors were: 5.6 x 10~ 4 for stage 7 (40 training pairs, 10 outputs per pair), and 7.9 x 10 -4 for stage 8 (50 training pairs, 20 
outputs per pair). 
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the relative size of the changes decreases, till by stage 
10 (the 5-qubit system's training to the full set of 56 
training pairs), changes in the parameters are just a few 
percent. Note that in the two figures, the relative sizes 
of the final numbers (root square percent change in K a 
is 3.6%; in Cab, 6.7%) are comparable. This is reminis- 
cent of the well-known neural network technique called 
"bootstrapping" [20] , in which knowledge of part of a pat- 
tern allows tentative inference about the rest; in a sense, 
information about a smaller system can, in part, general- 



ize to a larger. This makes us hopeful that this technique 
may be viable for use even on large systems. 

Figure [4] shows the average total entanglement (sum 
over all outputs) of the Wfc states as a function of k, 
for k = 1 to 5, as computed by the QNN using the 
trained 5-qubit system. Since the amounts of 3-way, 
4-way, and 5-way entanglements are zero on this scale, 
the total entanglement is equal to the sum of the pair- 
wise entanglements for each state. These numbers change 
only very slightly when we look at so-called "flipped" W 
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TABLE VII. Stage 9: five qubit QNN entanglement targets and calculated average outputs. Stage 9 added the five 4-way GHZ 
states to the training set of stage 8. RMS error = 1.4 x 10~ 3 (55 training pairs, 25 outputs per pair). 
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TABLE VIII. Stage 10: five-qubit QNN entanglement targets and calculated average outputs. Stage 10 added the one 5-way 
GHZ state to the training set of stage 9. RMS error = 2.4 x 10 -3 (56 training pairs, 26 outputs per pair). 
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states[in] instead (where, instead of a single excitation 
shared among N qubits, we have a single so shared - 
e.g., W 3 = -^(|110)±|101)±|011))). Recently Kimble et 

al.[2T] have succeeded in measuring the number of modes 
(k) sharing a single excitation photon; while the QNN is 
not, of course, calculating exactly this, our results do 
show similar separations among states of different k. 

Figure [5] shows the pairwise entanglements of the (nor- 
malized) states a|0001) + /3|0010) + |0100) + |1000), cal- 
culated by the QNN using the trained 4-qubit system, 
as a function of both a and f3. Results are not signifi- 
cantly different when using the trained 5-qubit system. 
When a = and /3 = 0, this is (|01) + |10>) <g) |00) = 
EPRab ® 1 00), which has full pairwise AB entangle- 
ment while all other entanglements are zero; our re- 



sults match this. When a = 1 and (3 = 0, this is 
1 0001) + 10100) + 1 1000) , which is the three-way W state in 
qubits ABD. In terms of pairwise entanglement, a three- 
way W state is equally pairwise entangled in all three pos- 
sible combinations AB, AD, and BD, so we would expect 
the red, green, and blue surfaces to come together at that 
point, while the yellow (AC), cyan (BC) and magenta 
(CD) surfaces are zero there. This is in fact what we see 
in the figure. Similarly, when a = and ,8 = 1, this is 
|0010) + |0100) + |1000), which is the three-way W state in 
qubits ABC; at this point the AB (red), AC (yellow), and 
BC (cyan) entanglements are all approximately equal and 
the AD (green), BD (blue), and CD (magenta) entangle- 
ments are zero. When both a — 1 and ,8 = 1, this is the 
four-way W state |001) + |0010) + |0100) + 11000), equally 
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TABLE IX. Final (stage 10) trained parameters for entan- 
glement of 5-qubit system, in MHz. Each row is a different 
parameter, as labelled. Each parameter is a function of time; 
the first row shows the timeslice number. 



t 


1 


2 


3 


4 


K A 


2.1733 


2.4021 


2.4688 


2.4139 


K b 


2.5842 


2.1591 


2.2702 


2.4710 


K c 


2.3984 


2.3252 


2.3501 


2.4294 


K D 


2.3893 


2.3128 


2.3451 


2.4327 


K e 


2.3926 


2.3144 


2.3446 


2.4274 


EA 


0.7443 


-0.1558 


-1.1248 


-0.7230 


EB 


0.9903 


-0.4483 


-1.0216 


-0.4294 


EC 


1.1027 


-0.4998 


-0.8678 


-0.4752 


EB 


0.7783 


-0.2453 


-1.1015 


-0.6125 


EE 


0.7703 


-0.2472 


-1.0941 


-0.6116 


Cab 


-0.3576 


0.2797 


-0.5974 


0.1521 


Cac 


-0.1761 


0.1800 


-0.5796 


-0.0604 


Cab 


-0.1687 


0.1802 


-0.5558 


-0.0471 


Cae 


-0.1563 


0.1758 


-0.5324 


-0.0402 


Cbc 


-0.2228 


0.1806 


-0.5462 


0.0340 


Cbb 


-0.2138 


0.1808 


-0.5227 


0.0431 


Cbe 


-0.2026 


0.1785 


-0.5014 


0.0484 


CCD 


-0.3052 


0.2375 


-0.4795 


0.1959 


CCE 


-0.2965 


0.2358 


-0.4638 


0.1801 


Cde 


-0.2896 


0.2265 


-0.4465 


0.1771 
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FIG. 2. Root square relative difference, AK^ / K, for 
both Ka and Kb, as a function of training stage. The sums 
are over the timeslices, and the differences are to each suc- 
cessive stage. Stage 1 is 2-qubit training]!?]; stage 2, 3-qubit 
pairwise; stage 3, 3q pairwise plus 3- way; stage 4, 4q pairwise; 
stage 5, 4q pairwise plus 3-way; stage 6, 4q pairwise plus 3- 
way plus 4-way; stage 7, 5q pairwise; stage 8, 5q pairwise plus 
3-way; stage 9, 5q pairwise plus 3-way plus 4-way; stage 10, 
5q pairwise plus 3-way plus 4-way plus 5-way. 
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FIG. 3. Root square relative difference in each of the param- 
eters ea, Eb and Cab as functions of training stage. Stages as 
m Figure [2] 



FIG. 1. Root mean squared error per training pair as a func- 
tion of epoch (pass through the training set), for the 5-qubit 
system on the pairwise training set (stage 7), summed over 
the 10 outputs per pair. 



entangled in all six pairwise ways AB (red) , AC (yellow) , 
AD (green), BC (cyan), BD (blue), and CD (magenta). 

The value for the pairwise entanglement of the three- 
way W states deserves some further comment. In the 
original training set for the two-qubit system, it was 
found [T7] that the value for the entanglement of the "P" 



state, 1 01} + 1 10) + 1 1 1 ) , that resulted in the smallest 
overall error (training plus testing), was 0.44: that is, us- 
ing this value, rather than some other, led to the great- 
est possible self-consistency for the QNN method. We 
therefore used this same value for P g) |0) for each of 
the three pairs, in training the three-qubit system [T§] . 
The three-way W state Wabc = |001> + |010) + |100) 
is not a "P" state in any of the three pairs, of course, 
but it is related in the following sense: we can write 
Wabc = |001) + (|01) + |10))® |0), that is, an EPR state 
in two of the qubits, plus a contamination term. In the 
same way, we can think of "P" as being a (fully entan- 
gled) EPR state, plus an amount of contamination, the 
inclusion of which diminishes the entanglement. Thus 
for self-consistency we would expect each of the pairwise 
entanglements of a W3 state to be approximately 0.44, 
as in fact occurs. 



flncangiemenc 01 w states 



Entanglement cf u|1 1C>H-|aD0>+ [i|111> 








FIG. 4. Total (pairwise) entanglement of Wfc, calculated by 
the QNN using the trained parameters listed in Table (IX} as 
a function of k. Each data point is the average of the total 

entanglement for the 2 fe_1 f ^ ^ states of the form Wfc ® 

|0) . Error bars show the standard deviation at each point, 
which increases from from 0.7%, for the Wb state, to 2.8% for 
the Wh state. The line is drawn to guide the eye. 
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FIG. 5. Pairwise entanglement of a| 0001) +/3|0010> + 1 0100) + 
1 1000), as a function of a and /3, as calculated by the QNN. 
Each color represents a different entanglement measure: red, 
the pairwise entanglement between A and B, Oab\ yellow, 
that between A and C, Oac, green, Oad\ cyan, ObC, blue, 
Obd; and magenta, Ocd- 



We can also look at states which have three-way 
(or higher) entanglement. Figure [6] shows the pairwise 
(red) and three-way (green) entanglement of the states 
a|110) + /3|111) + |000), calculated by the QNN using 
the trained 3-qubit system, as a function of a and /3, 
and compared with the residual or 3-tangle of Coffman|9 
(magenta). Results were not significantly different when 
using the trained 4-qubit or five-qubit systems. When 



FIG. 6. Entanglement of a|110) +/3| 111) + |000), as a function 
of a and /3, as calculated by the QNN, and compared with the 
3-tangle The red surface shows the pairwise entanglement 
between qubits A and B, Oab\ the green, three-way entangle- 
ment among A, B, and C,OabC, and magenta is the 3-tangle. 



a = and (3 — 1, this is the pure GHZ state for three 
qubits, and the three-way entanglement is maximal; both 
the 3-tangle and the QNN calculate this to be one. In- 
deed the QNN three-way entanglement (green) tracks the 
3-tangle (magenta) fairly well over the entire range of pa- 
rameters shown. This can also be seen in Figure [7J where 
we show the tracking of both pairwise and residual tangle 
by the QNN for states of the form |100)+/3|010)+7|00T). 

With the trained net, we can evaluate any pairwise or 
N-wise entanglement for mixed states, as well, without 
having to find an optimal decomposition]!]]]. Figure [8] 
shows the computed entanglement of states of the form 
a\GHZ 5 )(GHZ 5 \ + /3\W B )(W 6 \ + (1 - a - p)\W B )(\W B \, 
as a function of both a and f3. As expected from sym- 
metry, all the pairwise entanglements lie atop one an- 
other, so only one is shown, for clarity; all the three- 
way and four-way entanglements are zero. As this is 
a mixture, not a superposition, there is no interfer- 
ence among the three contributing states, and the \W B ) 
and \W B ) states contribute equivalently to the pairwise 
(red) entanglements. Compare this relatively feature- 
less behavior with that shown in Figure [9j which shows 
the calculated pairwise (red), four- way (blue), and five- 
way (green) entanglements of the superposition states 
a\GHZ 5 ) + P\W 5 ) + (1 - a - /3)\W 5 ). Again, all the 
pairwise entanglements (red) lie atop each other, by sym- 
metry, as do the three-way (not shown) and the four-way 
(blue). The three-way entanglements are still all zero; 
however, the four-way (blue) are not , because there are 
combinations of the coefficients such that the superpo- 
sition contains some nonzero amplitude of each of the 
\GHZ4) states; in neither | W5) nor \W$) are there any 
states with exactly three qubits in the same state, so 
there is no three-way entanglement. In both Figure [9] 
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Entanglement of |1C0>+[i|01 0+-,|001 > 
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FIG. 7. Entanglement of 1 100) + /3|010) + 7|001), as a func- 
tion of 7 and /3, as calculated by the QNN, and compared 
with analytical results. Each color corresponds to a different 
entanglement measure: red, pairwise AB entanglement, Oab', 
and blue, pairwise AC entanglement, Oac, both computed by 
the QNN. Compare these to tab, the tangle between A and 
B, in cyan; and tac, the tangle between A and C, in yellow 
(both computed analytically.) The residual tangle is identi- 
cally zero for all these states, in agreement with the QNN's 
calculated Oabc, in green. 



Entanglement mixture 
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FIG. 8. Entanglement of the mixed states a\GHZ 5 ){GHZ s \ + 
P\W 5 )(W s \ + (l-a-l3)\Ws}(\W5\, as a function of both a and 
p. Each color represents a different entanglement measure: 
red, the pairwise entanglement; blue, four-way entanglement; 
and green, five- way entanglement. 



and Figure [8] we see the expected limiting behavior: for 
a = 1 and f3 = 0, we have maximal five- way entanglement 
only; for a — and (3 — 1, we have maximal pairwise en- 
tanglement only, spread over all ten possible pairs (see 
also Figure |4j) 



FIG. 9. Entanglement of the superposition states a\GHZ^) + 
P\W 5 ) + (1 - a - P)\Ws), as a function of both a and /3. 
Each color represents a different entanglement measure: red, 
the pairwise entanglement; blue, four-way entanglement; and 
green, five-way entanglement. 

VI. CONCLUSIONS 

We have presented a new approach for the estima- 
tion of pairwise and multiqubit entanglement, applica- 
ble to any state, whether pure or mixed. Using dynamic 
quantum backpropagation learning, we have shown that 
a quantum system can be trained to compute its own 
degree of entanglement. No prior state reconstruction or 
tedious optimization procedure is necessary, nor is "close- 
ness" to any particular state. We envision an experi- 
mental implementation in which parameters could be, to 
begin with, roughly estimated from simulations such as 
those presented here (in the neural network literature, 
this is called "offline" or "batch" training [22].) These 
parameters could then be refined experimentally ("on- 
line" training of the neural network.) Because the state 
of the system at each timeslice must be known in order 
to do quantum backprop, it can only be done in simu- 
lation; some other parameter adjustment method would 
be needed experimentally, possibly a reinforcement or ge- 
netic algorithm. The experimental refinement procedure 
might also provide for the inclusion of features and in- 
teractions present in the physical state but not in the 
(simplified) model. In any case, once good values for the 
parameters are determined, any state's entanglement can 
be estimated by experimental measurement. And while 
we have only carried out simulations through a system 
size of 5 qubits, the additional training necessary signifi- 
cantly decreases as the size of the system grows, raising 
hopes for the method's feasibility even for large qubit 
networks. 

It should be emphasized, though, that this method pro- 
vides only an estimate of the entanglement (a "witness" , 
not a "measure.") Agreement with analytical results, 
while good, is not exact. Also, we have not explored here 
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application to phase shifted superpositions and mixtures, 
which present their own problems |10tjT2| IT7 ] . but which 
were not present in our training sets. Further work is 
needed in this area, and is ongoing. 

That being said, we believe our present work is still 
of great value, providing as it does a way of bypassing 
difficult analytical work in addressing the need[3T] for 
experimental measurements to determine entanglement. 
Moreover, our approach may have other important ad- 
vantages. Classically neural networks have proven fault 
tolerant and robust to noise; they are also famously used 
for noise reduction in signals. In quantum systems there 
is also the problem of decoherence. It is possible that 
quantum neural networks may be well suited for deal- 
ing with these types of problems in quantum computing. 



There are undoubtedly many other possible applications 
for the use of learning or AI methods as shortcuts to 
constructing quantum algorithms. 
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